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A precision measurement of the mass difference between the D° and D*° mesons has been made 
using 316 pb _1 of e + e _ annihilation data taken at yjs = 4170 MeV using the CLEO-c detector. 
We obtain AM = M(D*°) — M(D°) = 142.007 ± 0.015(stat) ± 0.014(syst) MeV, as the average 
for the two decays, D° — > K~n + and D° —> K~n + n~n + . The new measurement of AM leads 
to — 2006.850 ± 0.049 MeV, and the currently most precise measurement of the binding 

energy of the “exotic” meson X(3872) if interpreted as a D°D*° hadronic molecule, Aj,(X(3872)) = 
M(D°D*°) - M(X(3872)) = 3 ± 192 keV. 


Of all the claims and counterclaims for the so-called 
“exotic” mesons, which do not fit in the pictures of 
conventional qq mesons [l|, the most intriguing one is 
X(3872). Its existence has been confirmed from nu¬ 
merous measurements, by Belle [2|, CDF Q, DO 3, 
BaBar 0, LHCb 0 and CMS 0 , and its mass, width, 
and spin are respectively, M(X(3872)) = 3871.69 ± 
0.17 MeV, r(X(3872)) < 1.2 MeV, and J PC = 1++ 0. 
Although many different suggestions for the structure of 
X(3872) exist in the literature 00 , the closeness of 
the X(3872) mass to the sum of the masses of the D° 
and D*° mesons, and the smallness of its width have 
made the suggestion that it is a weakly bound hadronic 
molecule made of the D° and D*° mesons extremely at¬ 
tractive 0. To submit this provocative suggestion to 
experimental test it is important to measure the binding 
energy of X(3872), indeed to determine if it is bound at 
all. This Letter reports on the results of just such a mea¬ 
surement. Throughout this Letter we use the PDG Q 
convention for units, with masses in MeV and momenta 
in MeV/c, and inclusion of charge-conjugate states is im¬ 
plied. 

A measurement of the binding energy of X(3872) 
requires the knowledge of three masses, M(X(3872)), 
M(D °), and M(D*°), with the most accurately deter¬ 
mined value of AI(D*°) obtained by measuring the mass 
difference, AM = M(D*°) — M(D°). Since the dis¬ 
covery of X(3872) in 2003, the precision in the value 
of the mass of the X(3872) has steadily improved from 
±800 keV, originally, to the present average with error 
of ±170 keV [8j because of numerous improved measure¬ 
ments. Similarly, the precision of the value M(D°) has 
improved, from ±1000 keV, originally, to ±180 keV by 
a CLEO measurement of M(D°) in 2007 [l4} . and to 
±40 keV due to two recent hig her-precision measure¬ 
ments of M(Z?°) by BaBar [15[, and our recent pub¬ 
lication jTfij. As a consequence, the determination of 
the binding energy of X(3872) as a D°D*° molecule has 
changed from (600 ±600) keV in 2007 04!] to (126 ±204) 
keVin 2014 [IB]. 

Through all these improvements, the mass difference 


AM has remained fixed at the value, AM = 142.120 ± 
0.070 MeV as measured by CLEO in 1992 using data 
taken at the T(4S') resonance 17]. To determine the 
binding energy of X(3872) with the highest possible pre¬ 
cision, it has become imperative to make a new higher- 
precision measurement of AM = M(D*° ) — M(D°). 
In this Letter we report on such a measurement us¬ 
ing data taken at the ^(4160) resonance, which decays 
into D*W, D*D, and DD. We analyze D*° -T D°n°, 
and D° decays, D° —> K~n + (henceforth Ki r), and 
D° —> K~n + TT + TT~ (henceforth K3n). 

We use 316 pb" 1 of e + e _ annihilation data taken 
at = 4170 MeV with the CLEO-c detector. The 
CLEO-c detector 0 consists of a CsI(Tl) electromag¬ 
netic calorimeter, an inner vertex drift chamber, a central 
drift chamber, and a ring-imaging Cherenkov (RICH) 
detector, all inside a superconducting solenoid mag¬ 
net providing a 1.0 Tesla magnetic field. The accep¬ 
tance for charged and neutral particles is |cos0| < 
0.93. Charged-particle momentum resolution is cr p /p = 
0.6% @ 1 GeV/c. Photon energy resolution is (Je/E = 
2.2% @ 1 GeV, and 5% @ 100 MeV. The detector re¬ 
sponse was studied using a GEANT-based [IB] Monte 
Carlo simulation. 


We select events with well-measured tracks by requir¬ 
ing that they be fully contained in the barrel region 
(| cos 01 < 0.8) of the detector, and have transverse mo¬ 
menta > 120 MeV/c. 

In our previous article on the precision measurement 
of the mass of the D° meson |16(, we made a precision 
recalibration of the CLEO-c solenoid magnetic field and 
determined a correction of (2.9 ± 0.4) x 10~ 4 in the de¬ 
fault calibration of the CLEO-c magnetic field. The data 
we use in the present investigation were taken just after 
this recalibration. We use the same corrected field as 
determined in our previous article in the present Letter. 

Charged pions and kaons were identified using informa¬ 
tion from both the drift chamber dE/dx and the RICH 
detector. First, it was required that the dE/dx of the 
charged particle track be consistent within 3cr of the re¬ 
spective pion or kaon hypothesis. For tracks with mo- 
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FIG. 1. Reconstructed mass spectra for the candidates for 
the decays (a) D° —> Kir, (b) D° —> K3n. The verti¬ 
cal dashed lines show the regions in which we accept D° 
mesons, M(D° and ) = 1850 — 1880 MeV, and M(D° and ) = 
1855 — 1875 MeV, for D° —> K-k and D° —s> K3n channels, 
respectively. 




FIG. 2. Momentum spectra in data for (a) D° —> Kir and 
(b) D° — y K3tt for the candidates identified in Fig. 1. As de¬ 
scribed in the text, the enhancements correspond to momenta 
for the different expected final states with two charm mesons. 
The vertical dashed lines define the intervals of p(D° and ) in 
which the events were accepted for analysis. 


menta > 700 MeV, for which information from the RICH 
detector is available, the log-likelihood Lf ICH , as de¬ 
scribed in Ref. [lsj, was constructed. For such tracks, 
the combined dE/dx information a dE ^ dx , and the log- 
likelihood Lf ICH are used to distinguish between parti¬ 
cle hypotheses i and j , 


A L itj = [ a dE/dx f - (af /dx ) 2 + L f ICH - L f CH . (1) 


For tracks with momenta < 700 MeV, for which RICH 
information is not available, 

A Lij = (< j dE/d y - ( a dE/dx ) 2. (2) 

To identify pions, it was required that A< 0. For 
kaons, it was required that A ,K > 0. 

We reconstruct 7 r° —> 77 decays using photons only 
in the barrel region, |cos0| < 0.80. Photon candidates 


are defined as calorimeter showers with E 1 (barrel) > 30 
MeV and a transverse energy spread consistent with 
that of an electromagnetic shower. The photon candi¬ 
dates from 7 r° decays are required to have a two-photon 
invariant-mass, M{ 77 ), within ±15 MeV of the nomi¬ 
nal M(n 0 ) = 135.0 MeV Q, and to have neither pho¬ 
ton candidate combine with another photon candidate in 
the event to obtain an invariant mass closer to M(n 0 ). 
These candidates were kinematically fit with M ( 77 ) con¬ 
strained to the nominal 7r° mass in order to improve the 
energy resolution. 

We reconstruct single D*° candidates through the de¬ 
cays D*° 7 r°D°, D° Kit and D° -> K3ir, by first 
reconstructing the D° decay and then looking for the 
corresponding 7r° to make D*°. 

We select D° candidates using the standard CLEO D- 
tagging criteria, which impose very loose requirements 
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on the beam-energy-constrained D° mass, as described 
in Ref. [H| . 

In Fig. 1, we show the reconstructed invariant-mass 
spectra for D° —> Kn and D° —> K3n candidates. 

At -yds = 4170 MeV, D° mesons can be produced 
through several different decays: e + e“ —»• D*D*, D*D, 
and DD. As illustrated in Fig. 2, these different de¬ 
cays have different yields, and they populate different 
ranges of D° momenta. The enhancement in D° mo¬ 
menta at ft 500 MeV arises from e + e~ —> D*°D*°, and 
decays of both D*° and D*° to the corresponding D° 
and D°. The enhancement at ft: 720 MeV arises from 
e + e~ —> D*°D °, and D*° decaying to D°. The en¬ 
hancement at ft 950 MeV arises from e + e _ —> D°D°. 
GEANT-based [io| generic Monte Carlo is found to faith¬ 
fully reproduce both the observed enhancement, and the 
background, which arises from the other decays of D 
mesons. Events of interest for the present analysis are 
those in which at least one D*° (or D*°) is formed. In the 
momentum distributions of Fig. 2, there is an enhance¬ 
ment centered at p(-D° and ) ft 520 MeV, which arises from 
events with both D*° and D*°, which decay into Z?° and 
D °, respectively, and an enhancement at p(D® and ) ~ 720 


MeV which arises from events with only one D*° or D*°. 
The yield in the enhancement at p{D® and ) ft 520 MeV is, 
as expected, nearly a factor four larger than that in the 
enhancement at p(.D° and ) ft 720 MeV. 

We reconstruct D*° (or D*°) candidates for the events 
in the enhancements at p(-Dcand) ~ 520 MeV and 
p(Z)°and) ~ 720 MeV, and construct the distributions for 
AM = M(D* and ) — M(D° and ). These are shown sepa¬ 
rately for the decays D° —> Kn and D° —> K3n in Fig. 3. 

Since most of the instrumental uncertainties cancel in 
the measurement of the mass difference AM, the resolu¬ 
tion width of AM is approximately a factor of six smaller 
(FWHM ft 2.7 MeV) than that for the individual D° and 
D*° masses (FWHM ft 17 MeV). 

For both I\n and K3n, we fit the unbinned spectra 
to backgrounds parameterized by a second-order polyno¬ 
mial and signal parametrized as the sum of a Gaussian 
and another Gaussian with the same mean, but different 
widths on each side of the mean. Good fits are obtained. 

For clarity, we show the AM spectra for D° —> Kn 
(Fig. 3(a)) and D° —> I\3n (Fig. 3(b)), and the corre¬ 
sponding plots of residuals for the data binned in 250 keV 
bins. The results of the fits are listed in Table I. As is 
customary, the peak values from the unbinned fits are 
assumed to be measures of the corresponding mass dif¬ 
ferences, which are 
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FIG. 3. The AM spectra and fits for the decays (a) D° —> 
Kn, and (b) D° —> K3n. The fits (see text) were done to 
the unbinned data distributions, but for clarity the fits and 
corresponding residuals are shown for 250 keV bins. 


agree with those in Eqs. 3 and 4 within 20 — 50 keV, 
consistent with their statistical uncertainties. 


AM(I\n) =142.007 ± 0.018 MeV(stat), (3) 
AM(K3n) = 142.008 ± 0.027 MeV(stat). (4) 

We have analyzed the AM spectra for events in the 
enhancements at p(-D° and ) ~ 520 MeV, and p(D® and ) ft 
720 MeV, separately, and find that the results for AM 


We estimate systematic uncertainties in our results for 
AM from the following sources. The results are listed 
in Table II. The CLEO energy calibration for photons is 
based on the known n° mass, and on photon energies ex¬ 
tracted in the radiative decays i/>(2 S) —> 7 X c j(<7 = 1, 2), 
all of which are known with high precision. The uncer- 
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TABLE I. The results of the fits. 


Fit 

D u -4 Kn (Fig. 3(a)) 

D u -4 K 3tt (Fig. 

3(b)) 

Number of signal events 

6344 ±178 

3697 ± 134 


FWHM of AM distribution (MeV) 

ft 2.7 

« 2.6 


\ 2 /d.o.f of 250 keV binned fit 

1.10 

1.02 


AM from unbinned fit (MeV) 

142.007 ± 0.016(stat) 

142.008 ± 0.024 = 

7(stat) 

TABLE II. Systematic errors 

in AM = M(U c * a ° nd ) - 

- M(D c ° and ) in keV. 


Source 


D u -4 Kn D u 

-4 K'in 

Photon energy calibration (±0.4%) 


10 

10 


Charged particle momentum calibration (B-field ±0.4 x 10 -4 ) 
Signal shape (default vs. single Gaussian) 

Background shape (polynomial order +1) 

K ± mass (±16 keV) 

Cut variations in M(D)? and ), Fig. 1, (±2.5 MeV) 

Cut variations in p(Z)J? and ), Fig. 2, (±10 MeV) 

Cut variations in M(yy) (±2 MeV) 


3 

10 

1 

3 

4 
4 
4 


Total 


16 


16 


tainty in this calibration is estimated to be to ±0.4%. By 
rescaling the 7 r° photon energies by ±0.4%, we determine 
the resulting uncertainty in AM. In Ref. 16|, we made 
a precision determination of the CLEO solenoid B-field 
with an uncertainty of factor 0.4 x 10 -4 . This results 
in ±3 keV uncertainty in AM. The uncertainty in sig¬ 
nal shape is estimated using making alternate fits of the 
AM spectra by a single Gaussian in the restricted range 
of AM = 141 — 143 MeV. The systematic error in back¬ 
ground shape was obtained by increasing the order of the 
polynomial used in the fit by one unit. The PDG2014 
mass of the K± has an uncertainty of ±16 keV li¬ 
lt leads to ±3 keV uncertainty in AM. Contributions 
to the systematic uncertainties derived from variations 
of event selection requirements in M(D° and ) (Fig. 1), 
p(D® and) (Fig. 2), and M( 77 ) are seperately listed in 
Table II. Adding all systematic uncertainties in quadra¬ 
ture, the total systematic uncertainty is estimated to be 
±16 keV in AM(Kn) and AM{K3n). 

Our final results for AM, including systematic errors, 
are: 


AM(Kn) =142.007±0.018(stat)±0.016(syst) MeV,(5) 
AM(K3n) =142.008±0.027(stat)±0.016(syst) MeV.( 6 ) 

The results of the two D° decay final states are in good 
agreement (the agreement of the central values within 
1 keV is fortuitous). The weighted average of the two 
determinations, taking proper account of the correlations 
in the systematic errors in the two decay modes, is 

AM = 142.007 ± 0.015(stat) ± 0.014(syst) MeV. (7) 


The PDG2014 [s] average of AM based on the 1992 
measurement of Ref. 17 is AM = 142.12 ± 0.07 MeV 
based on a total of 1176 counts in Kn and K3n decays. 
Our measurement, based on 10,041 counts, has factor ft 
3.5 smaller overall uncertainty, and has 113 keV smaller 
value of AM. 

The present measurement (Eq. 7), AM = 142.007 ± 
0.021 MeV, and the latest average of D° mass, M(D°) = 
1864.843 ± 0.044 MeV, 0 lead to 

M(D*°) = 2006.850 ± 0.049 MeV. (8) 

This compares with the value M(D*% = 2006.96 ± 
0.10 MeV obtained by the PDG2014 @ from a simul¬ 
taneous fit of four masses, M(D* + ), M(D*°), M(D + ), 
and M(D°). 

Our measured masses lead to AI(D *°) + M(D°) = 
3871.693 ± 0.090 MeV. Using M(X{ 3872)) = 3871.69 ± 
0.12 MeV 0 we obtain the binding energy of X(3872) as 
a proposed D°D*° molecule, 

E b = (3871.693 ± 0.090) - (3871.69 ± 0.17) MeV = 

= 3 ± 192 keV. (9) 

The largest contribution to the uncertainty in the above 
result is due to ±170 keV uncertainty in the PDG2014 Q 
average value of the mass of X(3872). 

The negative limiting value of the binding energy E b 
implies that D°D*° system could be unbound by as much 
as 189 keV. The positive limiting value E b = 195 keV 
implies that the proposed D°D*° molecule, with reduced 
mass /i, has a minimum radius R = l/i/2 fiE b of 9.9 fm. 
Hopefully, our new result for the binding energy will shed 
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light on the continuing saga of the D°D*° molecule and 
other models of the structure of X(3872). 

This investigation was done using CLEO data, and as 
members of the former CLEO Collaboration we thank it 
for this privilege. This research was supported by the 
U.S. Department of Energy. 
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